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INTRODUCTION 


Further  improvement  in  chemical  rocket  technology  is  restricted  by  the  fundamental 
limitation  that  single-stage  devices,  independent  of  engine  efficiency,  will  always  have  negative 
payloads.  Therefore,  novel  concepts  that  might  serve  as  the  basis  of  new  propulsion  schemes  are 
of  great  interest.  One  concept  that  had  been  under  consideration  since  the  inception  of  the  HEDM 
program  is  that  of  metastable  molecular  fuels:  long-lived,  high-energy  content  molecular  species 
that  do  not  decay  by  radiation,  tunneling,  or  other  means  when  isolated  in  vacuum.  Consideration 
of  specific  impulse  virtually  requires  that  candidate  species  be  composed  of  first-row  atoms. 

The  goal  of  this  theoretical  research  program  is  to  evaluate  compounds  that  might  serve  as 
the  basis  of  new  propulsion  schemes.  The  program  was  originally  based  on  the  premise  that 
species  that  show  promise  as  high  energy  density  materials  may  exhibit  novel  bonding  mechanisms 
which  distinguish  them  from  conventional  stable  molecules. 

Two  possibly  overlapping  categories  of  molecules  were  defined  and  investigated:  (1) 
electron  deficient  compounds,  which  are  certain  compounds  of  B  and  Be  not  having  sufficient 
valence  electrons  to  distribute  two  per  chemical  bond  and  (2)  mixed  metal  clusters  of  the  form 

As  the  HEDM  program,  in  general,  and  this  research  program,  in  particular,  evolved, 
emphasis  has  shifted  from  seeking  new  metastable  compounds  to  serve  as  high  energy  fuels  by 
themselves  to  improving  on  conventional  liquid  hydrogen/liquid  oxygen  systems  by  devising  a 
way  to  incorporate  appropriate  additives  into  solid  hydrogen. 

The  boron  atom  has  been  identified  as  a  promising  additive  to  sH2.  In  studying  the  effects 
of  additives,  one  needs  knowledge  of  all  possible  compounds  of  boron  and  hydrogen.  Thus, 
calculations  of  other  compounds  of  boron  and  hydrogen  have  also  been  performed. 

The  classes  of  molecules  originally  proposed  included  compounds  based  on  "superalkalis” 
and  "superhalogens",  extraordinarily  stable  cations  and  anions  that  should  combine  to  give  species 
unusually  resistant  to  charge  transfer.  After  a  short  investigation,  it  was  determined  that  this 
category  was  unlikely  to  lead  to  compounds  of  sufficient  heat  of  formation  to  be  viable  HEDM 
candidates. 

This  report  describes  the  completed  study  of  the  electron  deficient  compound,  BH4,  as  well 
as  characterization  of  the  clusters  Li^B  and  Li^Be.  The  calculations  of  the  previously  unknown 
excited  states  of  BH2  undertaken  in  support  of  the  experimental  program  at  the  Phillips  Laboratory 
are  presented  and  a  brief  report  of  the  exploration  of  the  superalkali  concept  is  provided. 
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Theoretical  1^  calculations  tabulating  the  predicted  effects  of  including  these  compounds  as 
additives  to  the  liquid  hydrogen/liquid  oxygen  system  are  also  included. 
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STRUCTURE  AND  ENERGY  OF  THE  BH4  RADICAL 

Full  details  of  our  calculations1  of  BH4  are  given  in  a  paper  that  will  appear  in  the  Journal 
of  Physical  Chemistry  late  in  1993  and  is  attached  as  the  Appendix.  A  brief  summary  of  principal 
results  is  provided  here. 

In  this  work,  geometries  of  BH4,  the  dissociation  products  BH2  +  H2  and  BH3  +  H,  and 
the  possible  transition  state  for  the  addition  reaction  of  BH2  and  H2  were  optimized  at  the  MCSCF 
level  with  a  B(4s3p2d)H(4s2p)  basis  set.  As  illustrated  in  Table  1,  the  minimum  geometry  for  BH4 
is  predicted  to  have  C2V  symmetry  with  one  pair  of  short  B-H  bonds  and  one  pair  of  longer  B-H 
bonds.  The  structure  may  be  described  qualitatively  as  a  BH2  molecule  to  which  a  stretched  H2  has 
been  added  in  the  perpendicular  plane.  Analysis  of  the  wavefunction  indicates  that  of  the  seven 
valence  electrons,  four  are  contained  in  the  B-Hg  ana  B-Hb  bonds  while  there  are  only  three 
electrons  contained  in  the  longer  B-Hc,  B-Hd  bonds.  Comparison  with  calculations2  of  the  C2v 
symmetry  bidentate  form  of  the  ionic  compound  LiBH4  is  also  given  in  Table  1.  In  this  form,  the 
metal  ion  lies  in  the  plane.  The  principal  observation  is  that  the  longer  bonds  in  LiBH4 

are  shorter  than  in  BH4  due  to  the  additional  electron  in  the  HcBHd  moiety  as  compared  with  the 
radical. 


Table  1. 


BHaa 

LiBHa 

B-Ha 

1 .201 A 

1.200A 

B-Hc 

1.315A 

1.253A 

<Ha-B-Hb 

131.7° 

115.2° 

<Hc-B-Hd 

44.4° 

109.1° 

Hc-Hd 

0.994A 

2.041  A 

*MCSCF/TZP  (this  work)  bRamondo,  et  aL,  Ref.  2. 
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Harmonic  frequencies  for  BH4  and  LiBH4  are  listed  in  Table  2,  which  includes  a 
qualitative  description  of  the  modes.  Normal  modes  of  BH4  corresponding  to  the  BH2  symmetric 
and  asymmetric  stretches  (aj  and  t>2  symmetries  respectively)  can  be  easily  identified  Considering 
the  other  hydrogens,  Hg  and  Hd,  together  with  the  central  atom  as  a  BH2  unit,  the  corresponding 
modes  may  also  be  noted  with  reduced  frequencies,  as  expected  from  the  longer  bond  lengths. 


Table  2. 


CALCULATED  HARMONIC  FREQUENCIES 
(cm*1) 


BD4a 

BH4a 

LiBH4b 

ai 

1832 

2534 

BH2  sym  stretch 

2570 

1520 

2086 

BH2  (tong)  sym  stretch 

2133 

1003 

1413 

scissors  ( in  phase) 

1500 

751 

996 

scissors  (out  of  phase) 

1261 

bi 

1396 

1962 

BH2  (tong)  asym  stretch 

2059 

497 

634 

rock 

1170 

b2 

1999 

2656 

BH2  asym  stretch 

2612 

783 

1053 

twist 

1342 

a2 

628 

888 

twist 

1244 

BH2:  2456  sym  stretch 
2636  asym  stretch 
1 022  bend 

*MCSCF/TZP  (this  work). 
bRamondo,  at  al„  Ref.  2. 


Relative  energies  for  the  BH4  system  determined  from  multireference  single-  and  double¬ 
excitation  Cl  calculations  using  MCSCF  molecular  orbitals  are  reported  in  Table  3.  The  electronic 
energy  of  the  BH4  minimum  is  predicted  to  lie  19.23  kcal/mol  below  BH3  +  H  and  20. 10  kcal/mol 
below  BH2  +  H2.  Including  our  calculated  zero  point  energies,  the  BH4  species  is  bound  by 
14.83  and  14.55  kcal/mol,  respectively,  with  respect  to  BH3  +  H  and  BH2  +  H2.  Values  for  the 
heat  of  formation  of  BH4  and  dissociation  limits  at  0  K  and  298  K  are  also  listed  in  Table  3.  Wc 
have  chosen  to  base  these  numbers  at  0  K  on  a  value  of  22  kcal/mol  for  AHf°(BH3),  taken  from 
experimental3  and  theoretical4  work,  and  the  standard  value  of  51.6  kcal/mol  for  H.  Corrections  to 
298  K  were  made  following  the  procedure  given  by  Curtiss  and  Pople.4  Note  that  while  the  value 
in  the  JANAF  compilation5  for  AHf°(BH3)  of  26.4  kcal/mol  is  reasonably  similar  to  the  22 
kcal/mol  value  adopted  here,  the  JANAF  tables  give  an  incorrect  value  for  AHf°(BH2)  of  48.3 
kcal/mol  as  compared  with  73.3  kcal/mol  calculated  in  this  work. 
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Table  3. 


bh4  system  heats  of  formation 
(kcal/mol) 


Relative 

Energy/ZP 

AHf° 

AH2?8 

bh4 

-14.83 

58.5 

56.6 

bh3  +  h 

0.0 

73.6 

73.2 

BH3  +  H2 

-0.28 

73.3 

73.4 

bh2 

JANAF  entry3 

48.3 

aJANAF  thermochemical  Tables,  Third  Edition,  J.  Phys.  Chem.  Ref.  Data  14,  51  (1985). 
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EXCITED  STATES  OF  BH2 

As  stated  in  the  introduction,  the  boron  atom  has  been  identified  as  a  promising  additive  to 
sH2.  The  combustion  energy  of  B  is  second  only  to  that  of  Be,  among  the  first  row  atoms6  and  model 
calculations  by  Carrick7  have  predicted  addition  of  5  mole  percent  of  boron  atoms  to  sH2  would  result 
in  an  increase  in  theoretical  1^  of  80  seconds  or  21%  over  LH2/LO2. 

Accurate  predictions  of  the  1^  resulting  from  boron  atoms  or  small  molecules  containing  boron 
as  additives  to  solid  hydrogen  requires  accurate  heats  of  formation.  As  part  of  the  gas  phase 
spectroscopic  experimental  program  at  the  Phillips  Laboratory,  Carrick  and  Brazier  have  been 
investigating  the  spectrum  of  the  B2  molecule  produced  by  corona  discharge  of  diborane  in  helium  6,8 
It  was  originally  hypothesized  that  unidentified  lines  in  that  investigation  around  7000  A  could  be  due 
to  BH2,  although  it  was  later  shown8  that  they  can  be  attributed  to  B2.  While  the  two  lowest  states  of 
BH2,  the  X  2Aj  and  A  2Blf  have  been  characterized  experimentally9  and  theoretically10  we  were 
unable  to  find  any  previous  knowledge  of  the  higher-lying  excited  states  of  BH2.  Therefore,  we  have 
theoretically  predicted  the  excited  state  spectrum  of  the  BH2  molecule. 

The  results  of  our  investigation  are  presented  graphically  in  the  correlation  diagram  in 
Figure  1.  Relative  energies  with  respect  to  the  BH2  2Aj  ground  state  equilibrium  geometry  are 
given  in  Table  4.  In  constructing  this  diagram,  the  heat  of  formation  determined  in  our  study  of 
BH4  was  used  along  with  tabulated  atomic  levels  and  BH  molecular  excitation  energies  calculated 
in  this  work.  The  adiabatic  correlations  have  been  drawn  in  planar,  Cs  symmetry.  The  energetics 
give  the  BH2  ground  state  bound  by  83.6  kcal/mol  (3.62  eV)  with  respect  to  BH  +  H  including 
zero  point  energy  and  by  57.7  kcal/mol  (2.57  eV)  with  respect  to  BH  +  H. 

Table  4. 

RELATIVE  ENERGIES3  (eV)  OF  BH2  EXCITED  STATES 

Bent  Linear 

State  (R  =1.206  A,  8=127°)  State _ (R  =1.206  A) 

2Ai  0.0  ^aJ  0.386 

2Bi  1.262  2Bi  0.387 

2Ai,3S  5.724  2Ai,3s  4.671 

2B2  6.263  2Ai’  4.996 

2B2, 3p  6.612  2B2, 3p  5.604 

2Ai,3p  6.734  2A^3p  5.902 

2Bi,3p  6.836  2Bi'  3p  5.908 

2B2  7.359  2B2  6.645 

2Ai  7.559  2B2  7.878 

a)  With  respect  to  X2Ai  ,  equilibrium  geometry;  Energies  calculated  at  second  order  Cl  level 


ENERGY  (eV) 


BH2 

linear 


BH  +  H 


B  +  H2 


bh2 

bent 


BHI1!") 


BHf1!*  3s) 
BH(3I*  3s) 
BH(’a) 

BH(3!-) 

BH(’ri) 

BH(3nj 

BH(’I') 


CM  330M'  'S3 A 


Figure  1 .  Correlation  diagram  for  BH2.  States  are  connected  adiabatically  in  C$  symmetry  Bent 
geometry:  Calculated  equilibrium  geometry  of  the  X2Ai  state  (R«l  .206A,  0  -  127-; 
Linear  geometry:  R=1 .206A;  BH  asymptote:  R*1 .206A.  Energies  from  second 
order  Cl  calculations. 
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The  energies  displayed  here  were  determined  by  Second  Order  Cl  (SOCI)  calculations 
using  state-averaged  MCSCF  orbitals  with  the  basis  set  from  our  study  of  BH4  extended  by  two  3s 
and  two  3p  diffuse  functions  on  boron  and  one  diffuse  s  and  p  on  hydrogen.  In  each  molecular 
symmetry,  the  active  space  in  the  MCSCF  calculation  was  extended  by  the  appropriate  number  to 
accommodate  the  B  3s  and  3p  Rydberg  orbitals  of  that  symmetry.  The  energies  labeled  BH2  bent 
in  Figure  1  have  been  plotted  at  the  BH2  ground  state  geometry  (BH  bond  length  =  1.206  A  and 
bond  angle  =  127°).  At  this  geometry,  the  3s  Rydberg  state  lies  5.72  eV  above  the  ground  state  and 
the  three  components  of  the  3p  state  lie  between  6.61  and  6.84  eV.  However,  the  underlying 
BH2+  ion  has  a  linear  geometry.  One  would  therefore  expect  the  Rydberg  states  to  have  a  strong 
dependence  on  bond  angle.  This  effect  is  illustrated  in  Figure  2  where  the  excited  states  are  shown 
adiabatically  in  C2V  symmetry  as  a  function  of  bond  angle.  The  excitation  energies  at  the  linear 
geometry  (R(BH)  =  1.206  A,  0  =  180°)  are  also  shown  in  Figure  1.  All  of  the  Rydberg  states  are 
predicted  to  have  a  linear  equilibrium  geometry  as  are  the  valence  excited  states  with  the  exception 
of  the  third  2B2,  the  highest  state  plotted. 

We  have  also  studied  the  dissociation  path  (lengthening  one  BH  bond)  from  the  linear 
geometry  to  BH  +  H.  The  computational  procedure  included  state-averaged  CASSCF  calculations 
in  Cs  symmetry  (7  states  for  a'  symmetry  and  2  states  for  a"  symmetry)  followed  by  first  order  Cl 
calculations.  The  adiabatic  potential  curves  along  this  path  are  given  qualitatively  in  Figure  3. 
There  is  an  avoided  crossing  between  the  first  two  2A'  states.  Thus,  it  is  quite  possible  that  the  3s 
Rydberg  state  could  be  predissociated  to  BH  +  H.  The  accuracy  with  which  the  BH  excited  state 
energies  have  been  calculated  at  the  second  order  Cl  level  is  demonstrated  in  Table  5. 

Table  5. 

BH  ASYMPTOTES 


Calculated  Excitation  Energy  Te  (eV) 

BH  (eV)  Experiment/ 

State _ (R-1.206A) _ Literature 


x1!* 

0 

a3n 

1.30 

(1  21) 

Theory,  CISDa 

a 

2.95 

2.868 

Luh  and  Stwalleyb 

b3r 

4.70 

[4.66] 

Huber  and  Herzberg  b  ->  a 

C' 

5.90 

5.70 

Huber  and  Herzberg 

C  3I+  3S 

6.31 

[6.41] 

Huber  and  Herzberg  quote 
Pearson,  Bender,  Schaefer 

B  1l»3s 

6.42 

6.489 

Luh  and  Stwalieyb 

C  1L+ 

7.00 

6.86 

Huber  and  Herzberg 

aScuseria,  Geertsen,  and  Oddershede,  J.  Chem.  Phys.  90,  2338  (1989). 
bLuh  and  Stwalley,  J.  Mol.  Spectrosc.  102,  212  (1983). 
c  Pearson,  Bender,  Schaefer,  J.  Chem.  Phys.  55,  5235  (1971). 
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BH2  Energy  vs.  Angle  (FU1.206A) 


(A9)  A6jqu3 
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ENERGY  (eV) 


1  3  5  7  9  11 


RBH— H  (A) 

CM-330581  162A 


Figure  3.  Dissociation  of  linear  BH2  to  BH  +  H.  Potential  curves  are  drawn  adiabatically  in  Cs 
symmetry  as  a  function  of  the  BH  -  H  internuclear  distance.  Energies  from  first  order 
Cl  calculations.  2Ai  states  are  shown  by  solid  lines  and  2A”  states  are  given  by 
dashed  lines. 
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The  possible  transitions  of  the  BH2system  that  could  be  observed  experimentally  in 
emission  have  been  estimated  from  calculations  of  the  transition  moment  at  the  linear  geometry. 
The  crude  estimate  of  the  band  oscillator  strength,  f=2/3AE<r>2  was  used  along  with  the  standard 
relationship  between  oscillator  strength  and  Einstein  A  coefficient.  For  the  calculations  of  the 
transition  moment,  the  wavefunctions  were  recomputed  in  a  common  orbital  set.  The  allowed 
vertical  transitions  with  significant  strength  are  displayed  in  Figure  4  where  the  wavelength  in  nm 
and  Einstein  A  value  in  sec'1  are  shown.  Although  the  3s  Rydberg  state  is  likely  to  be 
predissociated,  the  rich  excited  state  structure  of  BH2  should  lead  to  emissions  that  may  be 
observed  experimentally. 

The  calculations  of  BH2  were  performed  in  collaboration  with  Dr.  Dahbia  Talbi  who  spent 
four  months  as  an  international  fellow  at  SRI  International  during  the  contract  period. 
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Vertical  transitions  in  emission  tor  linear  BH2(R=1 .206A,  0=180°).  Transitions 
labeled  by  wavelength  in  nm  (top  line)  and  A  value  in  sec'1  (lower  line),  estimated  as 
described  in  the  text. 
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MIXED  METAL  CLUSTERS:  Li3BE  AND  Li3B 


As  tabulated  in  the  following  section,  clusters  of  Li  with  beryllium  and  boron  are 
theoretically  predicted  to  improve  performance  when  used  as  additives  to  liquid  hydrogen. 
Calculations  were  performed  to  determine  the  structure  and  energy  of  the  Li3Be  system.  The 
computational  design  was  analogous  to  that  described  above  for  calculations  of  BH4  and  BH2. 
The  complete  active  space  for  geometry  optimization  included  5  electrons  distributed  in  4  orbitals. 
For  the  Q  calculation  of  energies,  the  active  space  was  enlarged  to  include  the  most  important  part 
of  the  2p  space,  as  determine  by  testing  of  the  full  2p  space. 

Optimized  geometries  are  reported  in  Table  6.  The  equilibrium  geometry  of  Li3Be  is  a 
planar  symmetric  species  in  agreement  with  the  previous  report  of  Pewestorf  et  al. 1 1 


Table  6. 

Li3Be  GEOMETRIES* 


Ll3Be  Minimum 

Ll'  2.7887  A 
Be 


Li' 


3.2038  A 


U2  Re«  2.9820 

Li  Be  (2Z+)  Re  -  2.5640 
a  Distances  in  A 


Li3Be 

Present  Work 

MC/DZ _ 

Li  Be  2.7787 

Li  Li’  3.1530 

Li’  Li’  2.5956 


Pewestorf  et  al.11 
HF/3-21G 
2.61 
3.13 
3.02 


U3  (2B2) 


Present  Work 
MC/DZ 

Li  Li’  2.8735 

Li’  Li'  3.2038 


Talbi  and  Saxon13 
MC/(5s2pld) 
2.820 
3.271 


Experiment  (Verma.  Koch,  Stwalley)14 
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Energies  for  the  Li3Be  system  are  listed  in  Table  7.  The  heat  of  formation, 

AHf  =  138  kcal/mol  reported  here  is  referenced  to  the  heat  of  formation  of  gaseous  atomic  Be 
from  the  JANAF  tables5  and  the  experimentally  derived  Li3  heat  of  formation  given  by  Wu. 12 

Table  7. 

RELATIVE  ENERGIES  WITH  RESPECT  TO  Ll3  +  Be  AND  HEATS 

of  formation  (kcal/mol)a 

LI3Be _ LI3  +  Be 

Relative  electronic  energy  -12.4  0.0 

Zero  point  energy  1.9  1.0 

AH?  137.9  149.5b 

a)  Assuming  AH? (Be)  *  76.5  kcal/mol  from  JANAF  tables,  Ref.5. 

AH?(Li)  -  73.0  kcal/mol  from  Wu,  Ref. 12. 

Geometry  optimization  of  Li3B  at  the  MCSCF  level  with  a  double  zeta  basis  set  gives  a 
diamond-shaped  symmetry  isomer  analogous  to  the  global  minimum  in  the  Li3H  system  and 
planar  D3h  symmetry  isomer,  analogous  to  the  BH3  geometry.  However,  the  lowest  energy 
isomer  is  a  T-shaped,  planar,  symmetry  species  that  lies  1.4  kcal/mol  below  the  symmetric 
D3h  structure  at  the  MCSCF  level.  These  isomers  have  also  been  investigated  with  the  MP2 
method.  Geometrical  parameters  arc  listed  in  Table  8. 

Table  8. 

LI3B  BONO  LENGTHS  (A) 


MC/DZ 

SCF/DZ 

MP2JDZ 

Djh  Symmetry 

B-Li 

2.204 

2.195 

2.233 

Li 

I 

#  imaginary 
frequencies 

(0) 

(0) 

(2) 

B 

/  \ 

U  Li 

MC/DZ 

MP2/DZ 

Symmetry 

B-Li 

2.419 

2.321 

U 

1 

B-Li’ 

2.282 

2.224 

ZLi’-B-li 

97.9 

92.4 

1 

_ B 

#  imaginary 
frequencies 

(0) 

(0) 

Li’  Li’ 


Li2  +  LI  Be 

8.5 

0.9 

157.9 


The  MP2  method  does  not  predict  stability  for  the  high  energy  diamond-shaped  isomer. 
Examination  of  the  MCSCF  wavefunction  indicates  sufficient  differences  from  the  single 
configuration  description  of  the  MP2  approach  to  lead  to  different  conclusions.  In  addition,  we 
determined  three  extrema  with  one  imaginary  frequency  that  lie  between  the  planar  diamond  and  the 
other  isomers  in  energy  at  the  MCSCF  level. 

Combining  these  results  with  our  previous  work  on  Li3H13  provides  an  interesting 
comparison.  Relative  energies  of  the  corresponding  structures  in  Li3H,  Li3Be,  and  Li3B  are 
displayed  graphically  in  Figure  5.  Cl  energies  are  given  for  Li3H9  and  Li3Be8  while  MCSCF 
results  are  used  for  Li3B.  Significant  contrasts  can  be  observed.  In  Li3B,  the  planar  diamond  is 
not  the  global  minimum  and  it  is  the  only  system  for  which  a  T-shaped  isomer  has  been  found. 

Earlier  in  the  HEDM  program,  interest  was  focused  on  high  energy,  metastable  isomers, 
but  when  use  of  atoms  as  additives  is  being  considered,  it  is  necessary  to  predict  the  lowest  energy 
compounds  that  could  be  formed  if  these  atoms  recombined.  It  is  clear  from  these  results  that 
various  stable  clusters  of  Li  are  to  be  expected  and  must  be  characterized  in  detail. 


15 


CM- 1060-1 

Figure  5.  Relative  energies  in  kcal/mol  for  U3H,  LisBe,  and  U3B  systems.  Energies  were 
determined  at  the  Cl  level  for  U3H  and  Li3Be  and  at  the  MCSCF  level  for  U3B,  as 
described  in  the  text. 
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ESTIMATED  THEORETICAL  I 


CALCULATIONS 


SP 

As  an  integral  part  of  program  planning,  theoretical  Isp  values  were  calculated  using  the 
AFAL  Theoretical  1^  Program  (micro  version)  for  various  compounds  of  boron,  lithium, 
hydrogen,  and  beryllium,  considered  as  additives  to  liquid  oxygen  and  liquid  hydrogen.  Results 
of  those  calculations,  carried  out  early  in  the  contract  period  and  modified  using  information 
obtained  later,  are  summarized  in  Table  9.  To  the  best  of  our  knowledge,  none  of  the  compounds 
listed  in  Table  9  have  been  characterized  experimentally.  Heats  of  formation  were  taken  from 
calculations  or  approximated  by  analogy  to  other  compounds.  The  value  for  B2H  provided  by  the 
calculations  of  Adams  and  Page15  should  be  quite  reliable.  Those  for  BH4  and  Li3Be  are  taken 
from  present  work.  The  heats  of  formation  of  Li3B  and  LiBg  should  be  viewed  as  guesses.  The 
L12B2  system  has  been  studied  by  Lammertsma.16  All  of  the  compounds  listed  in  Table  9  were 
theoretically  predicted  to  potentially  offer  an  advantage  over  the  LH2/LO2  system. 

Table  9. 


ESTIMATED  THEORETICAL  l8pa 


Additive 

AHf 

(kcal/mol) 

Additive 

Mole  Fraction 

O2  (liq) 

H2  (liq) 

lsp  (sec) 

none 

1 

4 

421 

B2H 

175 

1 

1 

10 

522 

4 

1 

15 

531 

3 

1 

16 

534 

2 

1 

16 

534 

bh4 

53 

2 

1 

9 

488 

2 

1 

8 

490 

2 

1 

7 

490 

U3B 

145 

1 

1 

16 

496 

ub6 

500 

1 

1 

16 

517 

1 

1 

18 

520 

400 

1 

1 

16 

487 

1 

1 

17 

490 

1 

1 

18 

492 

Li2B2 

162 

1 

1 

18 

489 

1 

1 

16 

489 

Li3Be 

138 

1 

1 

13 

513 

1 

1 

14 

512 

■  Calculated  with  AFAL  Theoretical  lsp  Program  (Micro  Version).  Density  of  additives  taken  as  1 .5  grams/cc 
Chamber  pressure  1000  psia.  Exhaust  pressure  14.696  psia. 
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"SUPERALKALI  -  SUPERHALOGEN"  COMPOUNDS 


On  the  basis  of  simple  molecular  symmetry  and  molecular  orbital  arguments,  Gutsev  and 
Boldyrev17  have  suggested  combinations  of  alkalis  and  halogens  or  chalcogens,  e.g.  Li2F,  Li30 
which  they  termed  "superalkalis"  because  they  are  expected  to  have  anomalously  low  ionization 
potentials  (IP's),  lower  than  that  of  the  isolated  alkali.  Similarly,  they  have  identified18 
"superhalogens",  e.g.  BeF3,  BF4,  with  anomalously  large  electron  affinities  (EA’s),  larger  than 
that  of  the  isolated  halogen.  Combinations  of  these  exceptionally  stable  cations  and  anions  might 
be  expected  to  lead  to  stable  ionic  species.  Charge-separated  species  have  been  of  interest  as 
HEDM  candidates  but  the  stability  of  the  ion-pair  species  investigated  previously  was  found  to  be 
limited  by  back-charge  transfer.  Thus  it  was  thought  these  unusually  stable  ions  could  lead  to 
interesting  possibilities  depending  on  the  heat  of  formation  of  the  resulting  compounds. 

The  first  step  was  to  ascertain  whether  the  qualitative  arguments17,18  based  on  approximate 
calculations  are  found  to  hold  with  more  rigorous  theoretical  methods.  Calculations  on  the 
molecules  B02  and  BF2  at  the  SCF  level  substantiated  the  arguments  based  on  orbital  energies. 

A  more  quantitative  study  was  undertaken  on  the  BF4  system,  which  was  selected  as  the 
first  superhalogen  candidate.  MCSCF  calculations  indicated  BF4‘  lies  more  than  100  kcal/mol 
below  BF3  +  F.  The  gas  phase  BF3  species  has  a  known5  heat  of  formation  of  -270.6  kcal/mol, 
and  that  of  the  F  atom5  is  given  as  18.4  kcal/mol.  These  values  indicate  that,  consistent  with  the 
"superhalogen"  argument,  BF4‘  is  an  extraordinarily  stable  anion.  However,  it  is  unlikely  any 
cation  could  be  found  with  a  heat  of  formation  high  enough  to  make  the  combined  charge- separated 
system  sufficiently  energetic.  In  addition,  the  recent  model  study  of  Fajardo19  points  out  possible 
limitations  on  the  densities  that  could  theoretically  be  achieved  for  charge  separated  species.  It  was, 
therefore,  determined  to  redirect  the  program  toward  the  studies  of  BH4,  BH2,  and  the  mixed  metal 
clusters  described  above. 
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THEORETICAL  INVESTIGATION  OF  THE  STRUCTURE  AND  ENERGY  OF 

THE  BH4  RADICAL 


Roberta  P.  Saxon 
Molecular  Physics  Laboratory 
SRI  International,  Menlo  Park,  CA  94025 

ABSTRACT 

Calculations  of  the  BH4  radical,  the  dissociation  products  BH3  +  H  and  BH2  +  H2,  and  the 
addition  reaction  of  BH2  and  H2  to  form  BH4  are  reported.  Structures  have  been  optimized  at  the 
MCSCF  level  with  a  polarized  basis  set.  The  Cj\  symmetry  equilibrium  geometry  may  be 
approximately  described  as  a  BH2  radical  with  a  stretched  H2  in  the  bisecting,  perpendicular  plane. 
Energies  have  been  obtained  by  large-scale  multireference  single-  and  double-excitadon  Cl 
calculations.  Including  zero-point  energy,  BH4  is  predicted  to  be  bound  by  14.55  kcal/mol  with 
respect  to  BH2  +  H2.  Heats  of  formation  at  0  K  and  298  K  are  also  reported.  Calculations  at  the 
Cl  level  suggest  the  addition  reaction  proceeds  without  an  energetic  barrier. 
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INTRODUCTION 


Electron-deficient  compounds  can  exhibit  a  variety  of  forms  due  to  the  rich  possibilities 
provided  by  large  numbers  of  empty  molecular  orbitals  and  thus  have  been  of  fundamental  interest 
in  investigations  of  chemical  bonding  mechanisms.  All  the  hydrides  of  boron,  for  example,  from 
BH  to  BH5  have  been  identified  directly  or  indirectly. 1-9  Boron  compounds  are  of  interest,  as  well, 
because  their  high  energy  content  makes  them  attractive  as  components  of  energetic  materials,10  yet 
the  actual  value  of  the  heat  of  formation  of  some  simple  boron  hydrides  is  still  a  matter  of  lively 
debate. 1'2*5-7  It  is  widely  recognized  that  the  standard  reference11  for  heats  of  formation  contains 
significant  inaccuracies  for  some  of  these  species. 

In  an  attempt  to  understand  the  energetics,  a  number  of  calculations  of  the  smaller  hydrides, 
BH,  BH2,  and  BH3  have  been  reported.1-5  At  the  other  end,  the  possible  existence  of  the  fully 
saturated  compound,  BH5,  proposed  as  an  intermediate  in  the  hydrolysis  of  BH4  has  attracted 
attention  for  a  number  of  years.8’9  Theoretical  studies  culminating  in  the  careful  work  of  Stanton, 
Lipscomb,  and  Bartlett9  have  established  that  BH5  is  an  intermolecular  complex  between  BH3  and 
molecular  hydrogen  and  that  the  very  modest  binding  is  due  entirely  to  electron  correlation  effects. 

The  open-shell  species  BH4  has  not  been  examined  with  the  same  intensity  although  an 
ESR  spectrum12*13  and  theoretical  geometry  predictions13-14  have  been  reported.  Election 
correlation  effects  noted  in  previous  studies  on  boron  hydrides2’9  may  be  expected  to  be,  perhaps, 
even  more  important  for  this  radical  species.  In  this  paper,  we  use  multiconfiguration  methods  to 
determine  the  structure  and  stability  of  the  BH4  radical.  Geometries  have  been  optimized  at  the 
MCSCF  level  and  the  energy  with  respect  to  dissociation  to  the  possible  dissociation  products, 

BH2  +  H2  and  BH3  +  H,  has  been  determined  by  MCSCF/CI  calculations.  The  reaction  path  for 
addition  of  H2  to  BH2  to  form  BH4  has  been  investigated  as  well. 
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CALCULATIONS 


Geometries  of  BH4,  the  dissociation  products  BH2  +  H2  and  BH3  +  H,  and  the  possible 
transition  state  for  the  addition  reaction  of  BH2  and  H2  were  optimized  at  the  MCSCF  level  using 
analytic  second  derivatives  at  each  step  to  facilitate  the  opdmization.  Calculations  for  the 
asymptotes  were  carried  out  on  the  supermolecule.  In  the  MCSCF  calculations,  seven  electrons 
were  distributed  without  restriction  in  the  complete  active  space  of  eight  orbitals  arising  from  B  2s 
and  2p  and  four  His  orbitals.  Geometry  optimizations  were  performed  first  with  a  double  zeta 
Gaussian  basis  set15  B(3s2p)H(2s)  and  refined  with  the  B(4s3p2d)H(4s2p)  basis16  used  for  the 
final  calculations  of  energies.  In  the  latter  basis,  the  scaled  triple  zeta  s  basis  of  Dunning17  for 
hydrogen  was  augmented  with  a  diffuse  s  function  with  exponent  0.09.  With  the  H(4s2p)  basis, 
the  SCF  energy  for  the  hydrogen  atom  is  0.499880  h  as  compared  with  0.492654  h  for  the  same 
basis  without  the  diffuse  function. 

Energies  were  determined  by  multireference  single-  and  double-excitation  Cl  calculations 
using  MCSCF  molecular  orbitals.  The  same  eight  orbital  active  space  as  for  the  MCSCF 
calculations  was  used;  the  lowest  occupied  orbital  corresponding  to  B  Is  was  frozen;  and  the 
remaining  orbitals  comprised  the  external  space.  The  configuration  list  included  all  single  and 
double  excitations  from  the  active  space  with  at  most  one  electron  in  the  external  space  and  all 
double  excitations  from  the  active  space  to  the  external  space  which  differed  by  no  more  than  two 
electrons  from  a  set  of  reference  configurations.  From  4  to  7  reference  configurations  were  used  in 
these  calculations  resulting  in  160,000  to  200,000  configurations  in  C2v  symmetry.  The  weight  of 
the  reference  configurations  in  the  Cl  wavefunctions  was  approximately  0.94  for  all  calculations. 
The  ALCHEMY  II  program  system18  was  used  throughout  the  study. 
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RESULTS  AND  DISCUSSION 


GEOMETRIES  AND  FREQUENCIES 

Geometrical  parameters  for  BH4  and  the  dissociation  products  BH3,  BH2,  and  H2 
optimized  at  the  MCSCF  level  are  given  in  Table  1.  The  identification  of  the  atoms  is  illustrated  in 
Figure  1.  The  global  minimum  on  the  potential  surface  for  BH4  is  predicted  to  have  C2v  symmetry 
with  one  pair  of  short  B-H  bonds  and  one  pair  of  longer  B-H  bonds.  The  structure  may  be 
described  qualitatively  as  a  3H2  molecule  to  which  a  stretched  H2  has  been  added  in  the 
perpendicular  plane.  The  bond  length  of  the  short  B-H  bond,  1.201  A  is  very  close  to  that  of  the 
isolated  BH2, 1.206  A,  and  the  Ha-B-Hb  bond  angle  is  also  similar,  131.7°  as  compared  with 
127.0°  for  BH2.  The  Hg-Hd  separation  of  0.994  A  may  be  compared  with  the  0.755  A  bond 
distance  in  H2  at  the  same  level  of  calculation. 

Harmonic  frequencies  for  BH4  and  dissociation  fragments  are  listed  in  Table  2,  which 
includes  a  qualitative  description  of  the  modes.  Normal  modes  of  BH4  corresponding  to  the  BH2 
symmetric  and  asymmetric  stretches  (ai  and  b2  symmetries  respectively)  can  be  easily  identified. 
Considering  the  other  hydrogens,  Hg  and  Hd,  together  with  the  central  atom  as  a  BH2  unit,  the 
corresponding  modes  may  also  be  noted  with  reduced  frequencies,  as  expected  from  the  longer 
bond  lengths.  (The  asymmetric  stretch  in  this  case  has  bi  symmetry  because  the  group  lies  in  a 
perpendicular  plane.)  All  other  modes  link  the  two  parts  of  the  molecule.  Harmonic  vibrational 
frequencies  for  BD4  are  also  reported  in  Table  2.  The  ratio  of  BD4  to  BH4  frequencies  varies  from 
0,71  to  0.78,  which  is  not  unexpected.  The  square  root  of  the  ratio  of  the  BD  to  BH  reduced  mass 
is  0.74. 

Geometries  and  frequencies  from  the  preliminary  determinations  with  the  double  zeta  basis 
set  and  results  of  previous  calculations  are  also  reported  in  Tables  1  and  2.  In  these  MCSCF 
calculations,  use  of  the  larger,  polarized  basis  set  had  only  a  modest  effect  on  the  predicted 


structures  and  frequencies.  The  greatest  differences  were  noted  for  the  B-Hc-Hd  subunit.  The  B- 
Hc  bond  length  was  reduced  by  0.013  A  and  the  Hc-B-Hd  bond  angle  was  reduced  by  7.1°  with 
the  final  basis  set  leading  to  a  reduction  in  the  Hc-Hd  separation  of  0.159  A.  Reasonable  agreement 
with  the  geometries  reported  previously13'14  for  BH4  and  for  the  smaller  hydrides  is  also  noted  in 
Table  1.  Previous  workers  did  not  list  harmonic  frequencies  for  BH4.  The  MCSCF  frequencies 
for  BH3  and  BH2  tend  to  be  as  much  as  10%  smaller  than  the  scaled  HF  frequencies1*5  which  is 
consistent  with  the  slighdy  longer  MCSCF  bond  lengths  resulting  from  the  modest  level  of  electron 
correlation  included  in  the  MCSCF  computational  model. 

The  MCSCF  wavefunction  for  BH4  is  strongly  dominated  by  the  single  SCF  configuration 
(Cl  coefficient  0.98).  Analysis  of  the  molecular  orbitals  indicates  that  of  the  seven  valence 
electrons,  four  are  contained  in  the  B-Ha  and  B-Hb  bonds  while  there  are  only  three  electrons 
contained  in  the  longer  B-Hc,  B-Hd  bonds.  The  singly  occupied  orbital  is  in  the  plane 
perpendicular  to  the  BH2  fragment  resulting  in  the  2Bi  symmetry  of  the  molecule.  Thus,  the  long 
B-H  bonds  in  BH4  are  real,  but  electron-deficient,  covalent  bonds. 

Results  for  BH4  determined  in  this  work  may  be  compared  with  values  predicted  for  the 
bidentatc  form  of  the  ionic  compounds  LiBH4  and  NaBH4.  The  equilibrium  geometry  of  the  BH4- 
negative  ion  is  tetrahedral  and  according  to  calculations19*20  and  experimental  observations21  the 
stable  configuration  of  UBH4  is  the  tridentate  form  in  which  the  BH4*  ion  is  only  modestly 
distorted  from  a  tetrahedron.  A  better  comparison  is  provided,  however,  by  the  symmetry 
bidentate  form  that  is  predicted19  to  be  a  transition  state  with  one  imaginary  frequency. 

Using  the  notation  of  Figure  1,  the  metal  ion  lies  in  the  Hc-B-Hd  plane  and  is  bonded  to 
both  He  and  H4.  Corresponding  bond  distances  and  bond  angles  are  tabulated  in  Table  3.  The 
principal  observation  is  that  in  UBH4  and  NaBH4,  the  longer  bonds  (labeled  B-Hc)  are  shorter 
than  in  BH4.  Assuming  these  compounds  are  really  ionic  in  nature,  i.e.,  Li+BH4\  there  is  an 
additional  electron  in  the  HcBH<j  moiety  as  compared  with  BH4  resulting  in  a  shorter  bond  length. 
This  also  results  in  a  larger  Hc-B-Hd  bond  angle  and  correspondingly  longer  Hc-Hd  bond  length  in 
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these  compounds.  Analogous  vibrational  frequencies  are  listed  in  Table  2.  They  may  be  observed 
to  be  somewhat  larger,  particularly  for  the  lowest  frequencies,  than  the  corresponding  frequencies 
in  BH4,  although  some  of  the  difference  may  be  accounted  for  by  use  of  different  theoretical 
methods. 

The  bonding  in  BH4  provides  an  interesting  contrast  to  the  tightly  bound  small  boron 
hydrides  at  one  extreme  and  to  the  situation  in  BH5  at  the  other,  where  there  is  no  bonding  between 
the  BH3  and  H2  fragments  at  the  SCF  level.  The  bond  length  of  1.315  A  in  BH4  is  intermediate 
between  the  ~1.2  A  B-H  bond  of  the  smaller  hydrides  and  the  1.47-1.53  A  B-H  bond  length,9 
depending  on  basis  set  and  type  of  calculation,  between  the  BH3  and  H2  segments  of  BH5.  Note 
also  that  the  increased  Hc-Hd  separation  of  0.994  A  in  BH4  also  supports  the  interpretation  that 
these  atoms  are  significandy  bound  to  the  central  atom  in  contrast  to  BH5  where  the  corresponding 
H-H  separation  is  0.75-0.79  A,  depending  on  the  calculation.9 

The  BH4  isomer  of  D2d  symmetry,  in  which  the  B-Ha  and  B-Hc  bond  lengths  and 
corresponding  bond  angles  are  equal,  was  also  investigated  with  the  double  zeta  basis  set  In 
agreement  with  the  previous  report,14  it  was  found  to  have  a  doubly  degenerate  imaginary 
frequency  and  to  lie  higher  in  energy  than  the  global  minimum  Cjv  structure.  It  is  therefore  not 
expected  to  be  chemically  significant  and  was  not  considered  further. 

ENERGIES 

Total  energies  for  BH4  and  the  dissociation  limits  as  well  as  relative  energies  with  respect 
to  BH3  +  H  are  listed  in  Table  4.  Comparing  results  at  the  MCSCF  level  for  the  preliminary 
geometry  optimized  with  the  double  zeta  basis  and  the  final  geometry  obtained  with  the  large  basis 
shows  a  very  modest  effect  on  the  energy  of  the  two  geometries,  corresponding  to  ~1  kcal/mol  for 
BH4.  Final  results  are  given  by  the  calculation  labeled  CID  which  includes  the  multireference 
Davidson  correction  for  the  effects  of  quadruple  excitations. 

The  electronic  energy  of  the  BH4  minimum  is  predicted  to  lie  19.23  kcal/mol  below 
BH3  +  H  and  20.10  kcal/mol  below  BH2  +  H2.  Including  our  calculated  zero  point  energies,  the 


BH4  species  is  bound  by  14.83  and  14.55  kcal/mol,  respectively,  with  respect  to  BH3  +  H  and 
BH2  +  H2.  These  results  may  be  compared  with  the  values  12.91  and  14.79  kcal/mol  reported 
previously14  from  MP4  calculations.  The  calculated  binding  energy  in  BH4,  although  modest  by 
the  standards  of  most  small  hydrides,  may  be  compared  with  the  largest  hydride,  BH5,  which  is 
predicted9  to  lie  5.4  kcal/mol  below  BH3  +  H2  in  electronic  energy,  which  reduces  to  0.9  kcal/mol 
when  zero  point  energy  is  taken  into  account.  This  is  consistent  with  the  shorter  B-H  bond  lengths 
in  BH4  than  in  BH5  as  noted  previously. 


Values  for  the  heat  of  formation  of  BH4  and  dissociation  limits  at  0  K  and  298  K  resulting 
from  the  energy  differences  in  Table  4  are  given  in  Table  5.  We  have  chosen  to  base  these  numbers 
at  0  K  on  a  value  of  22  kcal/mol  for  AH°(BH3)  and  the  standard  value  of  51.6  kcal/mol  for  H. 
While  there  remains  some  uncertainty  as  to  the  heat  of  formation  number  for  both  BH3  and  BH2, 
the  former  is  undoubtedly  more  reliably  known.  The  photoionization  studies  combined  with  careful 
analysis  of  several  different  reactions  of  Ruscic  et  al.7  provides  a  range  of  22.2  ±  3.4  to 
25.8  ±  1.7  kcal/mol  for  AH°(BH3)  while  theoretical  values  of  20.6  kcal/mol1  and 
22.3  ±  3.3  kcal/mol3  have  been  reported.  Corrections  to  298  K  were  made  following  the 
procedure  given  by  Pople  et  aL1  based  on  calculated  enthalpy  differences  for  BHn  and 
experimental  differences  for  the  atoms  B  and  H  in  the  standard  state.  This  results  in  AHf  (BH2) 
of  73.4  kcal/mol  which  may  be  compared  with  a  previous  theoretical  value1  of  74.8  kcal/mol.  As 
pointed  out  in  that  work1,  the  JANAF11  value  of  48  kcal/mol  is  based  on  old  appearance  potential 
data  subject  to  large  uncertainties.  Extrapolations11  are  based  on  an  assumed  linear  geometry, 
while  this  work  and  other  studies1*2  agree  on  a  bent  structure. 

The  heat  of  formation  of  BH4  at  0  K  and  298  K  is  predicted  to  be  58.5  and  56.5  kcal/moi 
respectively.  The  only  prev  ous  value  of  AHfCBHt)  in  the  literature,  to  our  knowledge,  is 


53.2  kcal/mol  calculated  within  the  MNDO  approximation.12 


ADDITION  REACTION 


Because  the  equilibrium  structure  of  the  BH4  molecule  has  a  C2v  geometry  that  resembles  a 
BH2  radical  with  a  slightly  stretched  H2  species  in  the  bisecting  plane,  the  addition  reaction  of  H2 
to  BH2  to  form  this  structure  is  of  interest.  Geometry  optimizadon  of  a  transition  state  for  the 
addition  reaction  analogous  to  those  for  BH4  and  the  dissociation  limits  was  performed  resulting  in 
the  structure  listed  in  Table  1.  At  this  geometry,  the  approaching  H2,  which  is  barely  stretched 
from  its  equilibrium  distance,  lies  in  the  plane  bisecting  the  BH2  radical.  The  distance  between  the 
closest  approaching  H  and  the  B  atom  is  1.844  A.  At  the  MCSCF  level,  the  transition  state  lies 
0.46  kcal/mol  above  the  BH2  +  H2  asymptote.  The  transition  state  is  characterized  by  an 
imaginary  frequency  of  448i  cm*1,  which  is  surprisingly  large  given  the  very  small  energy 
difference  between  the  transition  state  and  the  asymptote. 

When  the  energy  of  this  structure  was  determined  at  the  multireference  Cl  level,  however, 
the  electronic  energy  was  3.18  kcal/mol  below  that  of  the  asymptote.  The  "transition  state" 
structure,  located  at  the  MCSCF  level,  cannot  represent  the  location  of  a  barrier  on  the  true 
potential  energy  surface.  Even  when  the  zero-point  energy  of  this  point,  18.35  kcal/mol,  is  taken 
into  account,  the  "transition  state"  lies  0.4  kcal/mol  below  BH2  +  H2. 

In  order  to  see  whether  calculation  at  the  Cl  level  merely  alters  the  location  of  the  true 
transition  state,  we  also  performed  Cl  calculations  for  several  interpolated  points  in  between  the 
"transition  state"  and  the  asymptote  that  constitute  a  pseudo  reaction  path.  For  these  points,  the 
distance  between  the  approaching  H  and  the  B  was  fixed  and  the  other  coordinates  were 
interpolated  between  their  values  at  the  transition  state  and  the  asymptote,  defined  as  a  separation  of 
5A  between  the  approaching  H  in  H2  and  the  B  in  BH2.  The  energy  difference  between  an 
asymptotic  separation  of  5A  and  one  of  40A  is  0.025  kcal/mol  at  the  MCSCF  level.  Relative 
energies  at  the  MCSCF  and  Cl  levels  are  listed  in  Table  6.  It  appears  that  when  correlation  energy 
is  taken  into  account,  the  energy  of  the  system  goes  down  smoothly  as  the  hydrogen  molecule 


approaches  BH2.  These  calculations  suggest  that  there  is  no  energetic  barrier  to  the  addition  of  H2 
to  BH2. 

The  observation  of  a  small  barrier  in  a  radical  addition  reaction  "disappearing"  when 
reevaluated  with  a  larger  basis  set  and/or  degree  of  electronic  correlation  has  been  noted 
previously.  For  the  BH4  system,  Paddon-Row  and  Wong14  determined  a  similar  geometry  for  a 
transition  structure,  also  noted  in  Table  1,  which  was  found  to  be  slightly  lower  in  energy  than  the 
asymptote  when  recomputed  with  a  greater  degree  of  electronic  correlation.  Studies  of  insertion  of 
methylene  into  methane  and  silane^  and  of  halocarbene  cycloaddition^  also  reported  barriers  at 
the  SCF  but  not  at  correlated  levels  of  calculation. 


A-9 


SUMMARY  AND  CONCLUSIONS 


The  structure  and  energy  of  the  BH4  radical  has  been  studied  by  muiticonfigurarion 
theoretical  techniques.  The  equilibrium  geometry  of  BH4  is  found  to  have  C2v  symmetry  and  can 
be  approximately  described  as  a  BH2  radical  to  which  a  stretched  H2  has  been  added  in  the 
perpendicular  plane.  From  an  analysis  of  the  geometry,  frequencies,  and  MCSCF  wavefimction, 
all  of  the  bonds  in  BH4  may  be  considered  real,  covalent  bonds,  of  which  one  pair  is  electron- 
deficient  having  only  three  electrons.  This  is  in  contrast  to  the  well  studied8’9  BH5  system  which  is 
more  properly  described  as  a  complex  of  BH3  and  H2.  The  BH4  binding  energy  is  predicted  to  be 
14.55  kcal/mol  with  respect  to  BH2  +  H2  including  calculated  zero-point  energy  which  may  be 
contrasted  to  the  calculated9  binding  of  0.9  kcal/mol  of  BH5  with  respect  to  BH3  +  H2. 

The  addition  reaction  of  BH2  to  H2  to  form  the  C2v  equilibrium  geometry  of  BH4  was 
investigated.  Although  a  transition  stare  characterized  by  a  single  imaginary  frequency  and  a  very 
small  barrier  on  the  potential  energy  surface  was  determined  at  the  MCSCF  level,  subsequent 
investigation  at  the  Q  level  indicated  that  the  addition  reaction  takes  place  without  an  energetic 
barrier.  The  reduction  in  predicted  energy  barriers  when  treated  with  larger  basis  sets  and  greater 
inclusion  of  electronic  correlation  was  reported  in  other  studies,  as  well.  The  theoretical 
determination  of  energetic  barriers  in  radical  addition  reactions  appears  to  require  sophisticated 
calculations  for  reliable  results. 
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Table  1. 


OPTIMIZED  GEOMETRICAL  PARAMETERS 
(In  angstroms  and  degrees) 


BH4(C2v)2Bi 

MCSCFa 
(this  work) 

MP2/6-3lG(d,p)b 

UHF/6-31 1G**C 

B-Ha 

B-Hc 

ZHa-B-Hb 

ZHC-B-Hd 

Hg*Hd 

1 .201  (1 .207) 
1.315  (1.328) 
131.7  (128.7) 
44.4  (51.5) 

0.994  (1.153) 

1.117 

1.272 

129.2 

47.7 

1.029 

1.182 

1.288 

128.4 

48.4 

1.056 

BH2(C2v)2Ai 

MCSCF* 

(this  work) 

MP2/6-3lG*d 

HF/6-3lG(d)* 

Exptf 

B-H 

ZH-B-H 

1.206  (1.216) 
127.0  (128.4) 

1.188 

127.6 

1.185 

126.5 

1.181 

131 

MCSCF* 

(this  work) 

HF/6-3lG(d)* 

Expt9 

H-H 

0.755  (0.755) 

0.730 

0.742 

bh3  (DahjlAV 

MCSCF* 

(this  work) 

MP2/6-3lG*d 

HF/6-3lG(d)® 

B-H 

1.210  (1.217) 

1.191 

1.188 

Transition  State1* 
CsOA1) 

MCSCF* 

(this  work) 

MP2/6*31G(d,p)b 

B-Ha 

ZHa-B-Hb 

B-Hd 

Hc-Hd 

ZHc-B-Hd 

1.205  (1.211) 
128.5  (131.8) 
1.844  (1.689) 
0.766  (0.781) 
100.7  (103.4) 

1.180 

128.4 

2.069 

0.742 

97.1 

aResults  from  optimization  with  double  zeta  basis  in  parenthesis. 

bPaddon-Row  and  Wong,  Reference  14. 

cClaxton  et  al.,  Reference  13. 

dCurtiss  and  Pople,  Reference  2. 

ePople  et  al.,  Reference  1. 

fG.  Herzberg  and  J.W.C.  Johns,  Proc.  Roy.  Soc.  (London)  A298, 142  (1967). 

9K.P.  Huber  and  G.  Herzberg,  Molecular  Spectra  and  Molecular  Structure.  In  Constants  of  Diatomic 
Molecules  (Van  Nostrand,  New  York,  1979)  Vol  IV. 

hStructure  predicted  as  transition  state  at  MCSCF  level  (see  text).  Atoms  Hc  and  Hd  which  comprise  the 
approaching  H2  are  in  symmetry  plane  bisecting  Ha-B-Hb-  Hd  makes  the  closest  approach  to  B. 
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Table  2. 


CALCULATED  HARMONIC  FREQUENCIES 
(cm*1) 


BH4 

bh4 

BD4 

LiBH4b 

NaBH4b 

MCSCF 

MCSCF 

HF/6-31G** 

HF/6-31G** 

(this  work) 

(this  work) 

ai 

2534  (2535) 

BH2  sym.  stretch 

1832 

2570 

2533 

2086  (1941) 

BH2  (long)  sym.  stretch 

1520 

2133 

2153 

*1413  (1185) 

scissors  ( in  phase) 

1003 

1500 

1441 

*996  (611) 

scissors  (out  of  phase) 

751 

1261 

1278 

bi 

1962  (1770) 

BH2  (long)  asym  stretch 

1396 

2059 

2083 

*634  (742) 

rock 

497 

1170 

1164 

b2 

2656  (2670) 

BH2  asym  stretch 

1999 

2612 

2566 

*1053  (1043) 

twist 

783 

1342 

1286 

32 

*888  (781) 

twist 

628 

1244 

1240 

BH2 

MCSCF* 

(this  work) 

HF/6-3lG(d)c 

ai 

2456  (2460) 

sym.  stretch 

2728 

1022  (971) 

bend 

1128 

bi 

2636  (2631) 

asym  stretch 

2867 

h2 

MCSCF3 
(this  work) 

Expt^ 

4226  (4259)  4401 


bh3 

MCSCF3 
(this  work) 

HF/6-31  G(d)c) 

MBPT® 

Expt* 

i 

ai 

2455  (2436) 

sym  stretch 

2693 

2609 

aa 

1124  (1120) 

out-of-plane 

bend 

1225 

1185 

1140  88 

e' 

2579  (2587) 

stretch 

2813 

2756 

2808  (Ar) 

1201  (1147) 

1305 

1226 

’Modes  which  couple  the  two  parts  of  the  molecule. 
aResutts  from  optimization  with  double  zeta  basis  in  parenthesis. 
bC2y  symmetry,  bidentate  form,  Ramondo  et  al.,  Reference  19. 
cPople  et  at.,  Reference  1. 
dHuber  and  Herzberg,  Table  1 . 

eMBPT  with  B(3s2pld)  H  (2slp)  basis.  Stanton  et  al..  Reference  9. 
fM.  Jacox,  J.  Phys.  Chem.  Ref.  Data  19, 1387  (1990). 


A-1* 


Table  3. 

COMPARISON  OF  GEOMETRIES  FOR  BH4,  LIBH4,  and  NaBH4 
(in  Angstroms  and  Degrees) 


BH^a 

LiBHab 

NaBKib 

B-Ha 

1.201 

1.200 

1.204 

B-Hc 

1.315 

1.253 

1.249 

<Ha-B-Hb 

131.7 

115.2 

114.1 

<HC*B-H{j 

44.4 

109.1 

111.1 

Hc-Hd 

0.994 

2.041 

2.060 

a  MCSCF  (this  work). 

b  Ramondo,  et  al.,  Reference  19.  MP2,  6-31G*'. 
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Table  4. 


CALCULATED  TOTAL  ENERGIES  (hartrees)  ZERO-POINT  ENERGIES  (kcal/mol)  and 
RELATIVE  ENERGIES  (kcal/mot)  WITH  RESPECT  TO  BH3  +  H 

Total  Energies3 


MCSCF*5 

a 

CIDC 

Zero  Pt  Energy 

BH4 

-0.960307  (-0.958922) 

-1.054502 

-1.061982 

20.33 

bh3  +h 

-0.949652  (-0.949585) 

-1.026921 

-1.031334 

15.93 

BH2  +H2 

-0.945395  (-0.945284) 

-1.024962 

-1.029944 

14.78 

Relative  Energies*1 

MCSCF*5 

a 

CIDC 

BH4 

-6.69  (-5.86) 

-17.31 

-19.23 

BH2  +h2 

2.67  (2.70) 

1.23 

0.87 

aWith  respect  to  -26.0  hartrees. 

bFmal  geometry  (OZ  optimized  geometry). 

c  Cl  with  Davidson  correction  for  quadruple  excitations. 

d  Difference  in  electronic  energy  (without  zero-point  vibration). 
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Table  5. 


RELATIVE  ENERGIES  WITH  ZERO  POINT  WITH  RESPECT  TO  BH3  +  H 
AND  HEATS  OF  FORMATION  (kcal/mol) 

o  298 

Relative  energy/ZP  AHf  AH  f 


BH4 

This  work  -14.83 

MP2/6-31G(d,p)a  -12.91 

BH3+H  — 

BH2  +  H2 

This  work  -0.28 

MP2/6-31G(d,p)a  1-89 


58.5  56.6 

73.6  73.2 

73.3  73.4 


aPaddon-Row  and  Wong,  Reference  14. 


Table  6. 


RELATIVE  ENERGIES  (kcal/mol)  WITH  RESPECT  TO  BH2  +  H2 
FOR  ADDITION  REACTION 


MCSCF  a  CID0 


Transition  State”3 

B-Hd=»  1.84  A 

0.46 

-2.12 

-3.18 

B-Hd  »  2.9  Ab 

0.27 

-0.24 

-0.50 

B-Hd  «  3.9  Ab 

-0.02 

-0.10 

-0.12 

aSee  discussion  in  text. 

^Interpolation  procedure  for  geometries  described  in  text. 
CCI  with  Davidson  correction  for  quadruple  excitations. 
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FIGURE  CAPTION 


Figure  1. 


BH4  equilibrium  geometry.  (C2v  symmetry:  Ha  and  are  equivalent  and  Hc  and  H<j  are 
equivalent). 
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